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We present first results from three-dimensional radiation magnetohydrodynamic simulations of M-type dwarf stars with 
C05BOLD. The local models include the top of the convection zone, the photosphere, and the chromosphere. The results 
are illustrated for models with an effective temperature of 3240 K and a gravitational acceleration of log g = 4.5, which 
represent analogues of AD Leo. The models have different initial magnetic field strengths and field topologies. This first 
generation of models demonstrates that the atmospheres of M-dwarfs are highly dynamic and intermittent. Magnetic fields 
and propagating shock waves produce a complicated fine- structure, which is clearly visible in synthetic intensity maps in 
the core of the Ca II K spectral line and also at millimeter wavelengths. The dynamic small-scale pattern cannot be de- 
scribed by means of one-dimensional models, which has important implications for the construction of semi-empirical 
model atmospheres and thus for the interpretation of observations in general. Detailed three-dimensional numerical simu- 
lations are valuable in this respect. Furthermore, such models facilitate the analysis of small-scale processes, which cannot 
be observed on stars but nevertheless might be essential for understanding M-dwarf atmospheres and their activity. An ex- 
ample are so-called "magnetic tornadoes", which have recently been found on the Sun and are presented here in M-dwarf 
models for the first time. 
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1 Introduction 

Our picture of the atmospheric fine-structure of our Sun 
changed substantially during the past decades. Progress in 
instrumental performance resulted in unprecedented 
high-resolution observations that reveal a high degree of 
structure on a multitude of spatial and temporal scales. Such 
a complex phenomenon cannot be described in sufficient 
detail by means of otherwise elaborate semi-empirical one- 
dimensional model atmospheres (e.g. Vernazzaet al. 1981). 
The importance of small-scale spatial and temporal varia- 
tions in the solar chromosphere has been realized (cf . Ayres 
& Rabin 1996; Carlsson & Stein 1995; Solanki et al. 1991) 
and is now an essential feature of state-of-the-art numerical 
models. Furthermore, the layers of the Sun are no longer 
treated individually but as integral components of a dynam- 
ically coupled atmosphere (Wedemeyer-Bohm et al. 2009). 

Here we present a new class of 3-D models for M-dwarf 
stars that clearly demonstrates that their atmospheres are dy- 
namic and intermittent, too. This finding immediately gives 
rise to the question to what extent such inhomogeneities af- 
fect the emergent radiative intensity and related chromo- 
spheric diagnostics. In order to address this question, we 
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use model snapshots as input for spectrum synthesis calcu- 
lations and discuss the resulting implications of the spatial 
and temporal inhomogeneities for the interpretation of ob- 
servations and for the construction and usage of 
one-dimensional static model atmospheres. 



2 Numerical simulations 

The code. The numerical simulations are carried out with 
the 3-D radiation magnetohydrodynamic s code CO 5 BOLD 
(Freytag et al. 2012). The code numerically solves the equa- 
tions of (ideal) magnetohydrodynamics and radiative trans- 
fer together with a realistic equation of state and multi-group 
opacities. The simulation is advanced in time step by step 
and outputs a time sequence of 3-D snapshots and auxil- 
iary data. C0 5 BOLD is applied to a large range of different 
stars, among them the Sun (Steiner et al. 2008; Wedemeyer 
et al. 2004) and red giants (Freytag et al. 2002). Wende et 
al. (2009) used C0 5 BOLD to calculate local hydrodynamic 
models of M-type dwarf stars, which all include the top of 
the convection zone and the photosphere for a few selected 
parameter sets of T e ff and log g (see also Beeck et. al 201 1). 
The local simulations presented here include in addition a 
chromosphere and magnetic fields. The opacity tables are 
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a) Ca II K, AX = -3.93 nm 



Fig. 1 Visualisation of a 3-D model snapshot for a M- 
dwarf model with T e ff = 3240 K, log g = 4.5, and an initial 
magnetic field strength of \B\ = 100 G (2x2 checker- 
board). The bottom of the photosphere and the convection 
zone are grey-shaded, whereas regions of high magnetic 
field strength are coloured bluish. The red lines represent 
magnetic field lines. This image was produced with VAPOR 
(Clyne et al. 2007; Clyne & Rast 2005). 

based on PHOENIX NextGen data (Hauschildt et al. 1999). 
A particular challenge is that the typical computational time 
step for M-dwarf chromosphere models is on the order of 
only 1 ms, owing to the short time scales in the magnetized 
chromosphere. 

The models. A systematic model grid for different types 
of M-dwarf stars is currently under production. Here, we 
present exemplary models with T e ^ = 3240 K and log g = 
4.5, which represent (weak-field) analogues of AD Leo. The 
models have a horizontal extent of 1950 km x 1950 km and 
reach from -700 km in the convection zone to +1000 km in 
the chromosphere. All models start from a hydrodynamic 
simulation snapshot that has been relaxed from its initial 
conditions. The initial magnetic field, which is superimposed 
on the hydrodynamic model, is either unipolar and vertical 
or consists of regions of vertical field with opposite polar- 
ity (2x2 checkerboard). The individual models have initial 
magnetic field strengths ranging from 10 G to 500 G. 

3 Atmospheric structure and dynamics 

Under the conditions in the low photosphere outside strong 
magnetic flux concentrations (i.e., in quiet regions), the mag- 
netic field is effectively "frozen-in". Consequently, the flow 
field generated by the surface convection continuously re- 
arranges the magnetic field - like on the Sun. The initial 
magnetic field is swept into the intergranular lanes, where 
it is concentrated in form of continuous sheets and knots. 
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Fig. 2 Synthetic intensity images for a selected model 
snapshot. Top: CallK line wing (left) and line core (right). 
Bottom: continua at 0.3 mm (left) and 3.0mm (right). 



Especially in lane vertices, this process produces magnetic 
flux concentrations with field strengths in excess of 1 kG. 
The initial field is thus quickly transformed into a compli- 
cated non-uniform field. The complex geometry cannot suf- 
ficiently be described as a set of idealized "flux tubes". 

The magnetic field expands from these photospheric foot- 
points into the chromosphere, where it funnels out and fills 
the whole volume. The chromopheric magnetic field nev- 
ertheless exhibits a fine-structure that varies on short time 
scales. In the mixed polarity runs, magnetic loops form which 
connect footpoint areas of opposite polarity. The resulting 
3-D structure consists of entangled field lines which are of- 
ten twisted and rapidly change in time (see Fig.Q]). As a con- 
sequence, the model chromosphere varies strongly in space 
and time, resulting in strong variations in the physical state 
of the atmospheric plasma. Like for the Sun, a broad dis- 
tribution of gas temperatures is found in the upper layers. 
While upwards propagating shock fronts produce high tem- 
perature peaks in the chromosphere, the plasma in the post- 
shock regions cools to low temperatures. This non-uniform 
temperature distribution has to be taken into account when 
calculating meaningful average atmospheric stratifications. 



4 Synthetic intensity maps 

We illustrate the consequences of the complex atmospheric 
structure for the interpretation of observations. To this aim, 
we apply radiative transfer codes to three model snapshots 
for generating synthetic intensity maps, which can be inter- 
preted as virtual observations of the model. 
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Wavelength AX [nm] Wavelength AX [nm] 

Fig. 3 Synthetic intensity profiles of the Ca II K spectral 
line for a selected 3-D model snapshot. Left: The whole line 
profile. Right: Close-up of the line core. The grey area be- 
tween the 1 and the 99 percentiles marks the value range 
covered by the majority of the individual locations in the 
model. The corresponding 3-D average (thick black line, 
averaged horizontally over all columns in the 3-D model) 
exhibits a central emission peak. In contrast, the horizon- 
tally averaged model, i.e. the corresponding 1-D mean at- 
mosphere, results in a spectral line profile with no emission 
peak (red dot-dashed line). 

Calcium. First, we use MULTL3D (Carlsson 1986), which 
provides the detailed solution of the radiative transfer equa- 
tion in non-LTE (non-local thermodynamic equilibrium). The 
output contains intensity cubes (/ = f(x, A)) for 5 spec- 
tral lines (Ca II H, K, IRT) for each selected snapshot. Ex- 
emplary intensity maps in the line wing and in the line core 
are shown in Fig. [2 The line wing image essentially shows 
the granulation pattern originating from the low photosphere. 
The scene is completely different in the line core image. The 
interplay of magnetic fields and propagating shock waves 
causes an intermittent pattern of bright filamentary threads 
and dark regions, which mostly correspond to hot shock 
fronts and cool post-shock regions, respectively. The spec- 
tral line profiles accordingly vary significantly for the differ- 
ent locations. The intensity range covered by the individual 
Ca II K lines profiles is illustrated in Fig. [3 Despite the large 
spatial variations, the spatially averaged line profile still ex- 
hibits a strong central emission reversal peak. For compari- 
son, the model snapshot has been averaged horizontally, re- 
sulting in a one-dimensional model atmosphere. When com- 
puting the CallK line profile from this 1-D model atmo- 
sphere, again using MULTL3D, it does not show a central 
emission peak, which is in strong contrast to the full 3-D re- 




Fig. 4 Visualisation of a 3-D model snapshot for a M- 
dwarf model with T e ff = 3240 K, logg = 4.5, and an 
initial vertical magnetic field strength of \B\q = 100 G 
(2x2 checkerboard). The absolute magnetic field strength 
is color-coded in a plane that represents the bottom of the 
photosphere. The magnetic field lines (red) are concentrated 
in the intergranular lanes and funnel out in the atmosphere 
above. The plasma in the chromosphere follows the rotat- 
ing magnetic fields structures, resulting in spiral trajectories 
(green-blue lines). It is the first time that such a "magnetic 
tornado" has been found in a model of a cool star. 



suit. Similar 1-D average atmospheres for other model snap- 
shots can show emission in the line core but the differences 
to the 3-D result are always significant. It demonstrates that 
spatial averaging removes the strong fluctuations connected 
to the dynamic small- structure in the model chromosphere, 
resulting in synthetic spectral line profiles that carry only 
limited if not misleading information. 



(Sub-)millimeter continua. The continuum intensity at 
(sub-)millimeter wavelengths is another promising way to 
map the chromospheric layers. The large diagnostic poten- 
tial has been demonstrated for the Sun (Wedemeyer-Bohm 
et al. 2007). We use the radiative transfer code LINFOR3D 
(M. Steffen, http://www.aip.de/~mst/linj br3Djnain.html) to 
calculate the continuum intensity for wavelengths between 
0.3 mm and 9.0 mm, which lie in the range observable with 
the Atacama Large Millimeter Array (ALMA). Two exem- 
plary intensity maps are shown in Fig. [2 The average height 
in the atmosphere, where the intensity is formed, increases 
with increasing wavelength. For a plane-parallel stratified 
atmosphere, it would essentially map different chromospheric 
layers at the different wavelengths. However, the intermit- 
tent atmospheric structure results in large differences in the 
effective formation height range. At the same wavelength, 
it is possible to see deeper into the atmosphere at certain lo- 
cations than at others. The effective formation height range 
thus varies spatially and temporally and may make it diffi- 
cult to derive a meaningful atmospheric stratification from 
disk- averaged continuum intensity observations. 
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5 Implications for one-dimensional models 

Fuhrmeister et al. (2005) provide a semi-empirical model 
atmosphere for AD Leo, in which the gas temperature first 
decreases with height in the photosphere and increases es- 
sentially linearly in the chromosphere, followed by a char- 
acteristic jump in the transition region. The range of tem- 
peratures found in our 3-D model photosphere and with 
it the corresponding average matches the stratification by 
Fuhrmeister et al. (2005) quite well since the temperature 
fluctuations are only moderate in that layer. The temperature 
range drastically broadens in the 3-D model chromosphere, 
owing to the onset of shock formation. Peak temperatures 
on the order of 5500 K are found in many locations through- 
out the chromosphere, resulting in a height-dependent dis- 
tribution that overlaps with the semi-empirical model. The 
arithmetic horizontally averaged temperature, however, in- 
creases only mildly with height because hot and cold re- 
gions cancel. The temperature stratification obviously de- 
pends on how the temperatures enter the average. The rela- 
tion between (local) gas temperature and emergent intensity 
can be strongly non-linear for many spectral lines that are 
formed in the chromosphere, incl. the lines of Ca II. Deriv- 
ing a one-dimensional temperature stratification from spa- 
tially unresolved spectra can thus effectively lead to a higher 
weighting of the high temperatures connected to shock fronts 
Additional complications include largely varying for- 
mation height ranges, which makes it difficult to assign a 
definite height to an empirically derived temperature value. 
Consequently, the question arises how meaningful such 1-D 
average stratifications are. 

6 Small-scale dynamics 

Numerical simulations serve as a tool for analyzing and pre- 
dicting physical processes. The Sun plays an important role 
as a reference. The applied numerical methods can be tested 
by detailed comparisons of solar models with spatially re- 
solved observations of the Sun. For cool stars, adjustments 
and additional physics may be needed but nevertheless there 
are many qualitative similarities between solar and stellar 
models. A recent example are so-called magnetic tornadoes 
on the Sun (Wedemeyer-Bo?hm et al. 2012). They are gen- 
erated by the interplay of magnetic fields and vortex flows in 
the photosphere of the Sun. These vortex flows develop due 
to the conservation of angular momentum of cooled plasma 
in the downdrafts near the surface. In magnetically quiet 
Sun regions, the magnetic field is advected with the pho- 
tospheric flows, which forces the magnetic footpoints to ro- 
tate. The magnetic field mediates the rotation into the upper 
layers, where the situation reverses and the plasma follows 
the rotating field, resulting in spiral trajectories. The net en- 
ergy transport associated with these tornado-like events may 
significantly contribute to the heating of the solar corona. 

The models presented here suggest that photospheric 
vortex flows are abundant on the surface of cool M-dwarfs, 



too, like it was already found by Ludwig et al. (2002, 2006). 
As on the Sun, they force the magnetic footpoints to rotate 
and thus generate magnetic tornadoes in the chromosphere 
above. An example is illustrated in Fig.H It is plausible that 
these small-scale tornadoes are common on M-dwarf stars, 
too, which could be of importance for the energy balance of 
their outer atmospheric layers. 

7 Conclusions 

Already the first model generation clearly demonstrates that 
M-dwarf atmospheres are very dynamic and highly struc- 
tured on small spatial and temporal scales. Consequently, 
a detailed 3-D treatment is required for a realistic descrip- 
tion. One-dimensional model atmospheres, which are con- 
structed from a comparison to observed spectra, reproduce 
certain features but, by nature, cannot account for important 
properties of the atmospheres. Nevertheless, detailed 1-D 
models will remain complementary to 3-D models until all 
required physical processes can be treated numerically in 
sufficient detail. Such advanced 3-D models will be essen- 
tial for understanding the magnetic activity of M-type dwarf 
stars. 
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